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The quark-antiquark interaction potential is needed to compute properties of heavy quark 
systems, such as T mesons, or it production near threshold. In this paper, we will study 
the renormalization group improved potential at one-loop. The calculation makes use of 
NRQCD, formulated as an effective theory with an expansion in the velocity, v fllHl3 



The leading order term in v is the Coulomb potential, which has been computed to two- 
loop order 1)14,15] using QCD perturbation theory. The renormalization group running of 



this term is given by the QCD /3-function, as is well-known. We will compute the one- 
loop running of the order v and v 2 corrections to the quark potential, using a formulation 
of NRQCD introduced recently (T^j, and assuming mv 2 ^> Aq CD . In QCD the one-loop 
potential to order v 2 has been computed previously [fTTjnj]. For jj, — m the logarithms in 
the effective theory must agree with the logarithmic terms in these computations. We find 
agreement when some previously neglected terms are included in the spin-independent part 
of the quark potential. The renormalization group analysis allows one to resum logarithms 
of v in the effective theory, which is not possible using only the one-loop quark potential. 
The formalism we will use has been described in Ref. ||16|| , and will be called vNRQCD.[] 



In vNRQCD, one matches onto QCD at /i — m and then runs to lower scales in the effective 
theory using a velocity renormalization group (VRG). The VRG allows one to simultaneously 
sum logarithms of mv and mv 2 in the effective theory. In an alternative approach, called 
pNRQCD ||, the matching takes place in two stages, at fi = m and then at \i = mv. 
The logarithmic corrections to the potential were recently computed using pNRQCD by 



Brambilla et al. ||19|| . Our results agree with theirs when the resummed logarithms are 
expanded to linear order. 

The vNRQCD effective Lagrangian is written in terms of fields ip p which annihilate a 
quark, Xp which annihilate an antiquark, A£ which annihilate and create soft gluons, and A 11 
which annihilate and create ultrasoft gluons. The covariant derivative is = + igA^ = 
(D°, — D), so that D° = d° + igA°, D = V — igA, and involves only the ultrasoft gluon 
fields. The ultrasoft piece of the effective Lagrangian we need contains the quark, antiquark, 
and ultrasoft gluon kinetic energies, 




where the covariant derivative on ip p and Xp contain the color matrices T A and T A for 
the 3 and 3 representations, respectively. The coefficients in Eq. ([!]) do not run due to 
reparameterization invariance [poll . The potential interaction is 

Cp = -V a fSX T (p, p) 1p p >l 1p pf3 X-p'[ X-pr, (2) 

where a,/3, A,r denote color and spin indices. It is convenient to write the terms in V in 
matrix form. For example, 



1 Following a suggestion of A.H. Hoang. 
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FIG. 1. QCD diagrams for tree level matching. 
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represents the Coulomb interaction, and corresponds to Eq. (fj) with 
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The diagrams in Fig. |IJa,b give terms of the form 
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where k = p' — p and the notation of Titard and Yndurain [Tj| has been used: 
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The terms with coefficients V^, V^' 1 ^, V^ T \ and are order v 2 corrections to the lowest 
order Coulomb term, VjF\ Matching to the tree level diagram in Fig. |T]a at fi = m gives 

Vj T) = AiraJm) , Vi T) = Ana Am) , Vl T) = 



47ra s (m) 



vf* 
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(m) , vr»=-^#, v«=o. 



The annihilation diagram in Fig. |T]b gives the additional contributions 



Vg? = F 7ra.(m) 



(7) 



(8) 



These contributions have been separated from those in Eq. (^]) to facilitate comparison with 
results in the literature. In the color singlet channel there is no annihilation contribution. 

We have chosen to use the basis in which the potential V is written as a linear combination 
of 1 <8> 1 and T ®T. One can convert to the singlet and octet potential, using the linear 
transformation 



(9) 
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where C F = (iV 2 - 1)/(2N C ) and C A = N e . 

Potential terms can also be induced by operator mixing in the renormalization group 
flow. For instance, we will see that the running at one-loop induces order v 2 terms of the 
form 

y(T) y (l) 

V=(T A ®T A ) + \, (10) 

which are absent at tree level. Additional terms also occur when matching the potential at 
higher orders in the loop expansion. For instance, at one loop the order v correction 

V { k S) 7T 2 

kinglet = TTT (H) 

m |k| 

is generated. Ref. [18[ has = —C f (Ca — Cp/2) a s (m) 2 for the matching at /i = m. 

There are terms in the Lagrangian where ultrasoft gluons couple to potential operators. 
For instance, the diagram in Fig. |l]c generates the term 

C PU = J^yJ ABC (P " PO • (ffA°)[V T^ P X- P > j f B X- P ] . (12) 
The terms in the soft Lagrangian that we need for our computation are 

= E { \fK - <?xf + <p 9 i<p 9 + % q\ } (is) 
-g 2 E U V W, A «P$ V K>, ^ ^ 

P,P',<?,<?' 

+V [e*, q,]yW ^ P + (</y 1 T B Z^ Vp) (VVT^V)} + (V» - X, T - f ) . 



Here [/, W, Y", and Z are functions of (p, p', q, q') and matrices in spin and are generated by 
integrating out the intermediate offshell quarks and gluons in Fig. 0. The field c q is the soft 
ghost field, and (p q is the massless soft quark field with rif flavor components. The index a 
denotes the relative order in the v expansion. Performing the matching in Fig. ^| in Feynman 
gauge we find 

r(0) _ 1 rr(0) _ ( 2 P' " 2 P " <tf rr(0) _ (P " P' " ^ rr<P) 
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FIG. 2. The Compton scattering diagrams (a,b,c) generate the soft gluon coupling (d) in the 
effective theory. Diagrams analogous to (c) but with external ghosts or massless quarks generate 
the soft ghost and massless soft quark couplings in Eq. (13). 



u. 



(1) 



1.1 



u, 



(2) 



00 



U, 



(2) 



0; 



Tj(2) 



U. 



(2) 



Kl 



%cf e % i k cr k 
2m 



+ [25 ij q m + 5 <m (2p' - 2p - q) j + 5 jm (p - p' - qf] 



x 



(p+p') m + zc F e mfe '(p-p')V- 



2m(p' — p) 2 

c D (p' ~ P) 2 eg ia- ■ (p' x p) (p • q) 2 + (p' • q) 2 { (2 - c D )(p - p') • q 



i >m 2 q° 



Am 2 q° 



2m 2 (q°y 



Am 2 q° 



(1 ~ cp)q 2 
4m 2 g° 



[P • q (2p + q)' + P f -q (2p 7 - q) 1 ] ^ ic F [q x qf (p + p') • q 

4m 2 (g ) 2 

c D c 5 zcr ■ (p' x p) 



4m 2 (g ) 2 

| ( Cg -l)(p-p^ + q )» + (2p' - 2p - q) 



4m 2 



Am 2 



IP - P 



[p • q (p + p' + q) 1 + p' • q (p + p' - q) J ] *c F [(p - p' + q) x of ( P + p') ■ q 



4m 2 (g ) 2 
(eg - l)q* (p - p' - q) 



4m 2 (g )^ 



Am 2 



Am 2 



cd_ _c s mt ■ (p x p) 
2 



(P' - P) 2 



(P + pT (p + PT + 4(p-pQ-q^ + *c F ( P + p') J [(p - p') x a 



Am 2 q° Am 2 q° Am 2 q° 

ic F e ijk q k a- ■ (p + p') ic F e ijk cr k (p + p') ■ q (1 - cf^q^p - p' + q) j 



Am 2 q° 
(4 - c D )q 2 5^ 5 ij q° r 



4m 2 g° 



+ 



2m 2 



4m 2 g° 

cd _ eg ier ■ (p' x p) 
2 (p' - p) 2 



4m 2 q° 



5 



y(o) 

y(2) 
7 (o) 



0, 
1 
2m 



lp-pj-q 



(i) 



(P - P' + q) 



(i) 



-q 



-9 



2m(g ) 2 ' m 2mg° ' * u 2mg° 

q ■ (p + p') + ic F (T - [q x (p - p')] 



(i) 



2m 



: P ' - p) 2 ' 



y(i) 



2m(p' — p) 



c D g c 5 icr • (p' x p)g c 



5m' 



4m 2 (p' — p) 



r(0) 



:p' - p) 2 ' 



CD 

2 







4m 2 



c s zcr (p x p) 



IP " P 



0, Z 



A2) 



(1) 



-(P + p'Y-icfKp-p') x gf 
2m(p' — p) 2 



0. 



is not required for the calculation here. In Eq. (13) we have set p 2 = p' 2 , since this 
case is sufficient for our analysis. One could also set q° = |q|, however for calculational 
purposes we found it simpler to keep factors of q° explicit and ignore the energy poles this 
generates on the real axis. 

The coefficients of operators in the soft effective Lagrangian can run. In the soft regime 
the offshell quark propagators are identical to the quark propagators in heavy quark effective 



theory (HQET) |L0|,11]. Therefore, the divergence structure in the soft regime is the same 
as HQET. The running of C s can be computed indirectly using known results for the HQET 
Lagrangian, 

£ HQET = W^° + 5! + ^.B+ ^|[D.E] + ^f CT .(DxE-ExD)U. (15) 
I 2m 2m 8m z Sm z J 

The soft Lagrangian at fi = mv can be computed by first scaling the HQET Lagrangian to 
[i = mu, and then matching to the soft Lagrangian by computing the Compton scattering 
amplitude using HQET vertices in Fig. 0a,b,c. This gives the c F , Cd and c$ dependence in 
Eq. (|i~4"l). The running coefficients c F , cp, and cs were computed in Refs. |2T^p4|: 



c F [y) 
cd(v) 
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z -2C A /(3 + 



c s (v) = 2z~ Ca/i30 - 1 



(16) 



/20 32C F 
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Here z = a s (mz/)/a(m) and /5 = — 4T F nf/3, where T F = 1/2 is the index of the 

fermion representation. 

The renormalization group equation for the potential is 
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(17) 



where T is a 2 x 2 matrix and can be calculated as a power series in v and a s . The one-loop 
ultrasoft contributions to V are straightforward to compute. Assume that the potential has 
the form 
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v = (x A ® X A ) V (p, p') , 



(18) 



where X A is either T A or 1. An ultrasoft loop integrates over the ultrasoft loop- momentum, 
and due to the multipole expansion does not involve the labels p and p', so one can compute 
all one-loop divergent graphs with a single insertion of V. In Feynman gauge, the leading 
order graphs involve the ultrasoft vertex from the D° term in C u . The sum of all the 
graphs has no ultraviolet divergence, so there is no order v ultrasoft contribution to V [16 



Thus, the running of the quark potential involves mixing between different powers of v , i.e. 
running of the v 2 term proportional to the v term. To the order we are working, we need 
the ultrasoft vertices from the p • A/m operator, and insertions of p • V/m. Graphs with 
one insertion of V 2 / (2m) or p 4 / (8m 3 ) are the same order in the power counting but do not 
have ultraviolet log divergences. Graphs with one insertion of the operator in Eq. ( TT2] ) and 
one p • A/m vertex also do not contribute for this reason. The graphs which contribute to 
this mixing are listed in Table | and give 



p(us) 
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so Cd = (A c — 4/A c ). Here dimp is the dimension of the fermion representation. In the 
VRG method introduced in Ref. |L6||, the scale \i in ultrasoft loops is mu 2 , and the ultrasoft 
coupling constant is a s (mz/ 2 ). Therefore, the ultrasoft contribution to Eq. fllTD can be 
written as 



d_ 
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From Eqs. ( |T9"D and (|21| ) we see that the Coulomb potential induces running in V2 , V2 , 
and proportional to V^. 

In addition to the ultrasoft loops, one has the soft contribution shown in Fig. |3|. For 
the soft gluon loops all the infrared divergences are converted to ultraviolet divergences 
by tadpole diagrams, so with dimensional regularization all 1/e poles contribute to the 
anomalous dimensions. The divergent parts of the soft gluon, ghost, and quark loops in 
Fig. |3] give the running 



d_ 
^ dfi 



Bt®t 
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where for rif massless quarks we find 
14 
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TABLE I. Ultrasoft contributions to the running of the potential in Feynman gauge. The dot 
denotes an insertion of the potential in Eq. (|T8|), The p • A column gives the momentum factor 
from diagrams where the gluon coupling is due to the p • A interaction in the Lagrangian. The 
p • V column gives the factor from graphs in which the gluon vertices are from the D° interaction, 
and there are two insertions of the p • V operator on the quark lines in the loop. In d = 4 — 2e 
dimensions the ultraviolet divergent part of a diagram is — iV(p, p')a s /(2ire m 2 ) times the color 
and momentum factors. The sum of diagrams is gauge independent. 
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FIG. 3. Soft contributions to the running of the potential. The loop includes soft gluons, 

ghosts, and massless quarks. Graphs with two a = vertices from Eq. (O) gives the running of 

(T) 

the coefficient Vc ■ Graphs with one a = and one a = 1 vertex vanish, while those with one 
(7 = and one a = 2, or two a = 1 vertices contribute to the running of the order v 2 corrections 
to the potential. 



B 



6m 2 



C d + C A < 



22 (39 + 4c : 



3k 2 



6m 2 



19(p 2 + p' 2 ) 
3m 2 k 2 



(llc s + 10cp) 



r 2 



9m 2 ^ 1 
(c s + 2c F ) 



5C ^Tfk) 



+ 



T F n f f 1 (24-c D -i; 



2m 2 



A(p',p) 



r 2 

3m 2 



"\k 2 

r 2 
12m 2 



4m 2 



+ 



3m 2 
2m 2 k 2 



A(p',p) 



T(k) 



Here B 1(S1 and -Bt^t depend on the scale \i = mv through their dependence on cp, cd, 
and Cs- In the VRG, the soft coupling constant is a s {mv), so the soft contribution to the 
running is 



d 
t — 

dv 



Vt®t 



a 2 s {mv) 



B mx {mv) 
B T ®T{mv) 



(24) 



The total VRG running of the potential is given by adding Eqs. (^1|) and (|24]). The running 
of the Coulomb potential is proportional to the beta function, 



-2f3 a 2 s (mv). (25) 

Integrating Eq. (|25|) and using the boundary condition in Eq. (jfy gives the solution 

V c (u) = 4Tra s (mv) , (26) 

and large logarithms can be avoided by choosing v = \k\/m. This gives the expected 
result that the renormalization group improved Coulomb potential is given by choosing 
a s = a s (|k|). 

Order v corrections to the potential such as the one in Eq. (|ll|) do not run at one-loop. 
The non-trivial result is the renormalization group equations for the order v 2 terms in the 
potential. Combining Eqs. (|2~TD , (^4|), and using (|26|) we find 
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(27) 



Note that the soft contributions to the running depend on a 2 s {mv), and the ultrasoft contri- 
butions to the running depend on a s (mv)a s {mv 2 ). This is because the soft gluon coupling 
is g(mv), and the ultrasoft gluon coupling is g(mu 2 ). The ultrasoft gluon couples via a 
multipole-expanded interaction, so the ultrasoft interaction vertex does not involve momen- 
tum transfers of order mv. 

Our results can be checked by comparing with the one-loop formula for the color singlet 
quark potential in Ref. ||18|| . The nonrelativistic expansion of the QCD calculation has 
ln(|k|/m), ln(///m), and ln(A/m) terms, where a finite gluon mass A was introduced as an 



infrared regulator. By explicit computation of the box and crossed box diagrams |25| . |26 
found that the spin independent potential in Eq. (19) of Ref. 
termQ 



TBI] is missing the order v 



we 
2 



C fCa oi 



[P + P 

! k 2 
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(28) 



The ln(|k|) and ln(A) dependence in the full theory should be reproduced by the effective 
theory, so with /i — m all QCD logs should be reproduced. We find that the log(A/m) 
terms are reproduced by the ultrasoft diagrams in Table [I] and that the log(|k|/m) terms 
are reproduced by the soft loops in Fig. |3|. 

Solving the VRG equations in Eq. ( P?D with the tree level boundary conditions in Eq. ([7[) 
givesQ 



2 Note that in the calculation in Ref. [17| an expansion was made in (p + p') 2 /k 2 , so that the term 
in Eq. ( p8| ) was dropped. 

3 In NRQED the results are simpler since the coupling in the effective theory does not run. 
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(29) 



where 2 = a s (mu)/a(m). The annihilation contributions can be accounted for by adding 
the expressions in Eq. (Rp to Eq. (25p. 

The logarithmic corrections to the coefficients in the quark potential were considered 



by Brambilla et al. |T9| using pNRQCD, but were not resummed. Since only the leading 
perturbative logarithms were included, we can compare our results to theirs by expanding 
the resummed logarithms in Eq. (p9|). Brambilla et al. give expressions for the V's in the 
color singlet channel containing terms of the form 



a s (r), a s a s (r)ln(mr), a s a s (r) ln(/xr) . 



(30) 



For the soft contributions letting a s (r) — > a s (rrw) and ln(r) — > \a{l/mv) in Ref. |B| gives 
agreement^ with expanding the soft contributions in Eq. ( p9|) about z = 1 (recalling that 
z — 1 — fio/ (2ir)a s {mv) \n{mv/m)). For the ultrasoft contributions letting ln(/i/m) — > 
ln(mz/ 2 /m) in Ref. |19| gives agreement with the ultrasoft contributions in Eq. (|29[ ) with the 
expansion 
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(31) 



It would be interesting to compare our resummed coefficients to the analogous results com- 
puted in pNRQCD. Recall that there is an important distinction between pNRQCD and 
vNRQCD. In pNRQCD the scales mv and mv 2 are treated independently, while vNRQCD 
builds in the fact that the scales mv and mv 2 are not independent. Thus, in vNRQCD, the 



4 When comparing the expansion of the color singlet combination v| — CfV% the replacements 
Cd — > 8Cf — 3Ca and C\ — > CfCa/2 — Cp are also necessary to obtain agreement. 
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scaling of soft terms from m to mv, and ultrasoft terms from m to mv 2 , occurs simultane- 
ously, and \n(mv) and \n(mv 2 ) terms are resummed together. 

To see the effect of the running on the value of the coefficients in the potential, consider 
the case of top quark production near threshold. Projecting onto the color singlet channel 
gives the singlet coefficients = — CpV^. Using a s (m t ) = 0.108, the tree level 
values in Eq. (0) are: 

V r (s) = -1.81 , Vi s} = , V { s s) = 0.60 , Vi s) = 0.15 , V t W = 2.71 . (32) 

For a Coulombic system the velocity v is determined by solving a s (mv) = v. Using m t = 
175 GeV and the one-loop running of a s (fi) with rif = 5 gives v = 0.145. For v = v the 
running coefficients in Eq. (|29|) are: 

Vi s) = -1.49 , V^ s) = 0.63 , Vi s) = 0.53 , v| s) = 0.16 , V { t s) = 3.11 . (33) 

The most substantial change is to the value of which was zero at tree level. Using the 
results of Brambilla et al. with mr = 1/v and /x = mv 2 gives 

VW = -1.78 , V ( 2 s) = 0.68 , V { s s) = 0.53 , V| s) = 0.16 , V t (s) = 3.15 , (34) 

indicating that resummation of the logarithms has the largest effect on V^. For the re- 
maining coefficients the resummed values in Eq. (|33"D are fairly close to the coefficients with 
perturbative logarithms in Eq. (j34|). 

We would like to thank M. Luke, A. Pineda, and I. Rothstein for helpful discussions. 
This work was supported in part by the Department of Energy under grant DOE-FG03- 
97ER40546. 
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